CLAIMS 



1. A system estimation method for making state estimation 
robust and optimizing a forgetting factor p simultaneously in an 
5 estimation algorithm, in which 

for a state space model expressed by following expressions: 
x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
10 here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 
15 z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

a maximum energy gain to a filter error from a disturbance 

weighted with the forgetting factor p as an evaluation criterion 
20 is suppressed to be smaller than a term corresponding to a 

previously given upper limit value y f/ and 

the system estimation method comprises: 

a step at which a processing section inputs the upper limit 
value y f , the observation signal y k as an input of a filter and 

25 a value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
with the upper limit value y f ; 

30 a step at which the processing section reads out an initial 

value or a value including the observation matrix H k at a time from 
the storage section and uses the forgetting factor p to execute 
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a hyper Hoo filter expressed by a following expression: 

x x k |k = Fk-ix'Vilk-i + Kg, k (yk - HkFk-iX^k-xjk-i) 
here, 

x\ )k ; an estimated value of a state x k at a time k using observation 
5 signals y 0 to y k/ 

K s ,k; a filter gain, 

a step at which the processing section stores an obtained 
value relating to the hyper H*, filter into the storage section; 

a step at which the processing section calculates an 
10 existence condition based on the upper limit value y f and the 
forgetting factor p by the obtained observation matrix Hi or the 
observation matrix Hi and the filter gain K S/i/ and 

a step at which the processing section sets the upper limit 
value to be small within a range where the existence condition 
15 is satisfied at each time and stores the value into the storage 
section, by decreasing the upper limit value y f and repeating the 
step of executing the hyper filter. 

2. The system estimation method according to claim 1, 
20 wherein the processing section calculates the existence condition 

in accordance with a following expression: 

S'l = S-l l + 1 -^HjH i >0, i = 0,...,* (17) 

3. The system estimation method according to claim 1, 
25 wherein the processing section calculates the existence condition 

in accordance with a following expression: 
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-qE { + P7/>0, z = (),••, A 



here, 



i ^2 - pH{K si 



(18) 
(19) 
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4. The system estimation method according to claim 1, 
wherein the forgetting factor p and the upper limit value y f have 
a following relation: 

0 < p = 1 - x(Yf) ^ !' where x(Yf) denotes a monotonically damping 
function of y f to satisfy =1 and %(<») = 0. 

5. The system estimation method according to claim 1, 
wherein at the step of executing the hyper H«> filter, 

the processing section obtains the filter gain K s<k by 
following expressions: 



x k \ k = F k -ix k -i\k-i + K s<k {y k - H k F k-ix k -i\ k -i) 
K 3>k = Sk^iHliHkSk^Hl + p)- 1 
Zk\k = 2 k \k-i - Zk\k-\C T k R~l,CkZk\k-\ 
S M \k = (F k Zk\kFl) IP 



(10) 

(11) 
(12) 

(13) 



here, 



e fti = Zi\i- HiXu &o|o = £o, i7i| 0 -^o 

p 0 



R ek = R k + Ck2 k \k-iC T k , R k 
0 < p = 1 - x(7/) < 1, 7/ > 1 

-Fk^k\kF\ 



. 0 -pyj 



C k = 



(14) 
(15) 



Q kQ T = X(7/) P r rT 



(16) 
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wherein a right side of the expression (16) can be more 
generalized, 
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here , 

x k : the state vector or simply the state, 
y k : the observation signal, 
z k : the output signal, 
5 F k : the dynamics of the system, 
H k : the observation matrix, 

x^ k | k : the estimated value of the state x k at the time k using the 
observation signals y 0 to y k , 

X^kik: corresponding to a covariance matrix of an error of x^ k | k , 
10 K S/k : the filter gain, 

e f<j : the filter error, and 
R e , k : an auxiliary variable. 

6. The system estimation method according to claim 5, 
15 wherein the step of executing the hyper Hoc filter includes: 

a step at which the processing section calculates the filter 
gain K S/k by using the expression (12) based on an initial condition; 

a step at which the processing section updates a filter 
equation of the Hoo filter of the expression (11) ; 
20 a step at which the processing section calculates Z^k|k and 

Z^ k+ i|k by using the expression (13); and 

a step at which the processing section repeatedly executes 
the respective steps while advancing the time k. 



25 7. The system estimation method according to claim 1, 

wherein at the step of executing the hyper Hoo filter, 

the processing section calculates the filter gain K s , k by 
using a gain matrix K k and by following expressions: 
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*k\k = «*-i|fc-i + K sJc (y k - Hkik-nk^) (20) 
tf SiJfc = #*(:> l)/i? e ,*(l,l) , iiC* = P^p-^fc/Zri Jf 1 ) Jifljfc (21) 
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Where, 

R k = RlJiRl Rl = 
Re t k = Rk + Ck£k\k-\Cl, C k - 



' 0 




' 1 0 








0 ^27/ 


.0-1 



, Re,k = Rl k J l Rl k , x 0 |o = io (23) 



0(k) denotes a J-unitary matrix, that is, satisfies 0(k)J0H(k) T 
= j / j = ( J x ® I ) , I denotes a unit matrix, K k ( : , 1) denotes a column 
vector of a first column of the matrix K k , 

wherein Ji" 1 and Ji can be deleted in the expressions (21) 
and (22), 
here, 

x~ k | k : the estimated value of the state x k at the time k using the 

observation signals yo to y k , 

y k : the observation signal, 

F k : the dynamics of the system, 

K S/k : the filter gain, 

H k : the observation matrix, 

X^lk*. corresponding to a covariance matrix of an error of x^ k)k , 
0(k): the J-unitary matrix, and 
R e , k : an auxiliary variable. 



8. The system estimation method according to claim 7, 
wherein the step of executing the hyper H«> filter includes: 

a step at which the processing section calculates K k and 
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Z^k + i[k 1/2 by using the expression (22) ; 

a step at which the processing section calculates the filter 
gain K s , k based on the initial condition and by using the expression 
(21) ; 

a step at which the processing section updates a filter 
equation of the Hoo filter of the expression (20) ; and 

a step at which the processing section repeatedly executes 
the respective steps while advancing the time k. 
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9. The system estimation method according to claim 1, 
wherein at the step of executing the hyper Hoo filter, 

the processing section obtains the filter gain K s , k by using 
a gain matrix K k and by following expressions: 

&k\k = £jk-i|fc-i + KsAVk - H k x k -i\k-i) (61) 
K Stk = K k (:, 1) . Kk = pHK k RJ)Rl k (62) 



0 



0 



0 



C k +\L k R r jt 
P^L k R~l 



Q{k) (63) 



here, 0(k) denotes an arbitrary J-unitary matrix, and C k = Cjt+i^F 

is established, 

where 



R k — R£ J i R£ , R£ — 



?3 0 

0 ps^f 



Jl = 



1 0 
0 -1 



1 . I 

5 y 2 



Re,k = Rk + Ck£k\k-lC{, Cjfc = 



H k 



(23) 



here, 

20 x^ k | k : the estimated value of the state x k at the time k using the 
observation signals y 0 to y k , 
y k : the observation signal, 
K s , k : the filter gain, 
H k : the observation matrix, 
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0(k): the J-unitary matrix, and 
R e ,k: an auxiliary variable. 

10. The system estimation method according to claim 9, 
5 wherein the step of executing the hyper H«> filter includes: 

a step at which the processing section calculates K~ k by using 
the expression (63); 

a step at which the processing section calculates the filter 
gain K~ s , k based on the initial condition and by using the expression 
10 (62) ; 

a step at which the processing section updates a filter 
equation of the H* filter of the expression (61) ; and 

a step at which the processing section repeatedly executes 
the respective steps while advancing the time k. 

15 

11. The system estimation method according to claim 1, 
wherein at the step of executing the hyper filter, 

the processing section obtains the filter gain K s , k by using 
a gain matrix K" k and by following expressions: 
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K,, k = p*K k (:,l)/R e ,k(hl) 

-p-*L k R;lL k c k+i 



Kk+i 




0 


0 







0 



L J 

i*r,*+l = #r,* - L k C kn R ek C M L k 

Where , 



(25) 
(26) 

(27) 

(28) 

(29) 
(30) 



, = «(* + 1 - iv)] = K* + 1) u*], = Ki) ) o,... 1 o] 
P o 



1 0 
0 0 
0 1 



-1 0 
0 



, ^ 1 ,o = diag{p 2 J ^ ) ...,^ +2 } ) p = \- X {u) 
. ^o = 0, x 0 |o = xo, Kk=p~ h K k (31) 



wherein the above expressions can be arranged also with 
respect to K k instead of K~ k , 
here, 

y k : the observation signal, 

F k : the dynamics of the system, 

H k : the observation matrix, 

x A k |k: the estimated value of the state x k at the time k using the 
observation signals yo to y k , 

K s , k : the filter gain, obtained from the gain matrix K~ k , and 
R e/k , L~ k : an auxiliary variable. 

12. The system estimation method according to claim 11, 
wherein the step of executing the hyper H« filter includes: 

a step at which the processing section recursively 
calculates K" k +i based on a previously determined initial condition 
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and by using the expression (27) ; 

a step at which the processing section calculates thesystem 
gain K s , k by using the expression (26) ; 

a step at which the processing section calculates the 
5 existence condition; 

a step at which the processing section updates a filter 
equation of the Hoc filter of the expression (25) when the existence 
condition is satisfied, and repeatedly executes the respective 
steps repeatedly while advancing the time k; and 
10 a step of increasing the upper limit value y f when the 

existence condition is not satisfied. 



13. The system estimation method according to claim 1, 
wherein an estimated value z v k | k of the output signal is obtained 

15 from the state estimated value x\ )k at the time k by a following 
expression: 

z\|k = H k X* k | k . 

14. The system estimation method according to claim 1, 
20 wherein the Hoo filter equation is applied to obtain the state 

estimated value x^ k | k , 

a pseudo-echo is estimated by a following expression: 

7V-1 

dk = Y.hA k W-u k = 0,1,2,.-. (34) 

and 

25 an echo canceller is realized by canceling an actual echo 

by the obtained pseudo-echo. 



15. A system estimation program for causing a computer to 
make state estimation robust and to optimize a forgetting factor 
30 p simultaneously in an estimation algorithm, in which 
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for a state space model expressed by following expressions: 
x k+1 = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
5 here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 
10 z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

a maximum energy gain to a filter error from a disturbance 

weighted with the forgetting factor p as an evaluation criterion 
15 is suppressed to be smaller than a term corresponding to a 

previously given upper limit value y f , and 

the system estimation program causes the computer to 

execute : 

a step at which a processing section inputs the upper limit 
20 value y f , the observation signal y k as an input of a filter and 
a value including an observation matrix H k from a storage section 
or an input section; 

a step at which the processing section determines the 
forgetting factor p relevant to the state space model in accordance 
25 with the upper limit value y f ; 

a step at which the processing section reads out an initial 
value or a value including the observation matrix H k at a time from 
the storage section and uses the forgetting factor p to execute 
a hyper Hoo filter expressed by a following expression: 
30 x^ k | k = Fk-ix'Viik-i + K s , k (y k - H k F k -iX Vi|k-i) 
here , 

x^ k | k ; an estimated value of a state x k at a time k using observation 
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signals y 0 to y k , 

F k : dynamics of the system, and 
K S/ k; a filter gain, 

a step at which the processing section stores an obtained 
5 value relating to the hyper Hoc filter into the storage section; 

a step at which the processing section calculates an 
existence condition based on the upper limit value y f and the 
forgetting factor p by the obtained observation matrix H ± or the 
observation matrix Hi and the filter gain K s ,i; and 
10 a step at which the processing section sets the upper limit 

value to be small within a range where the existence condition 
is satisfied at each time and stores the value into the storage 
section by decreasing the upper limit value y f and repeating the 
step of executing the hyper filter. 

15 

16 . A computer readable recording medium recording a system 
estimation program for causing a computer to make state estimation 
robust and to optimize a forgetting factor p simultaneously in 
an estimation algorithm, in which 
20 for a state space model expressed by following expressions: 

x k+i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
here , 

25 x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 

z k : an output signal, 
30 F k : dynamics of a system, and 

G k : a drive matrix, 

a maximum energy gain to a filter error from a disturbance 
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weighted with the forgetting factor p as an evaluation criterion 
is suppressed to be smaller than a term corresponding to a 

previously given upper limit value y f , and 

the computer readable recording medium recording the system 
5 estimation program causes the computer to execute: 

a step at which a processing section inputs the upper limit 
value y f/ the observation signal y k as an input of a filter and 
a value including an observation matrix H k from a storage section 
or an input section; 
10 a step at which the processing section determines the 

forgetting factor p relevant to the state space model in accordance 
with the upper limit value y f ; 

a step at which the processing section reads out an initial 
value or a value including the observation matrix H k at a time from 
15 the storage section and uses the forgetting factor p to execute 
a hyper Hoo filter expressed by a following expression: 
x^ k | k = F k -iX%-i| k -i + K s , k (y k - H k F k _iX^ k _i| k _i) 
here , 

x^ k j k ; an estimated value of a state x k at a time k using observation 

20 signals y 0 to y k/ 

F k : dynamics of the system, and 
K g , k ; a filter gain, 

a step at which the processing section stores an obtained 
value relating to the hyper Hoo filter into the storage section; 

25 a step at which the processing section calculates an 

existence condition based on the upper limit value y f and the 
forgetting factor p by the obtained observation matrix H± or the 
observation matrix Hi and the filter gain K 3ti ; and 

a step at which the processing section sets the upper limit 

30 value to be small within a range where the existence condition 
is satisfied at each time and stores the value into the storage 
section by decreasing the upper limit value y f and repeating the 
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step of executing the hyper H*, filter. 



17. A system estimation device for making state estimation 
robust and optimizing a forgetting factor p simultaneously in an 
5 estimation algorithm, in which 

for a state space model expressed by following expressions: 
x k+ i = F k x k + G k w k 
y k = H k x k + v k 
z k = H k x k 
10 here, 

x k : a state vector or simply a state, 

w k : a system noise, 

v k : an observation noise, 

y k : an observation signal, 
15 z k : an output signal, 

F k : dynamics of a system, and 

G k : a drive matrix, 

a maximum energy gain to a filter error from a disturbance 

weighted with the forgetting factor p as an evaluation criterion 
20 is suppressed to be smaller than a term corresponding to a 

previously given upper limit value Yf, 

the system estimation device comprises: 

a processing section to execute the estimation algorithm; 

and 

25 a storage section to which reading and/or writing is 

performed by the processing section and which stores respective 
observed values, set values, and estimated values relevant to the 
state space model, 

the processing section inputs the upper limit value y f , the 
30 observation signal y k as an input of a filter and a value including 
an observation matrix H k from a storage section or an input section, 
the processing section determines the forgetting factor p 
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relevant to the state space model in accordance with the upper 
limit value Yf , 

the processing section reads out an initial value or a value 
including the observation matrix H k at a time from the storage 
5 section and uses the forgetting factor p to execute a hyper H«> filter 
expressed by a following expression: 
x"jc| k = Fk-iX^k-ijk-i + K s , k (vk - H k F k -iX^k-i|k-i) 
here , 

x^k|k; an estimated value of a state x k at a time k using observation 
10 signals y 0 to y k , 

F k : dynamics of the system, and 

K s ,k; a filter gain, 

the processing section stores an obtained value relating 

to the hyper Hoo filter into the storage section, 
15 the processing section calculates an existence condition 

based on the upper limit value y f and the forgetting factor p by 

the obtained observation matrix Hi or the observation matrix H ± 

and the filter gain K s<i , and 

the processing part sets the upper limit value to be small 
20 within a range where the existence condition is satisfied at each 

time and stores the value into the storage section by decreasing 

the upper limit value Yf repeating the step of executing the 

hyper filter. 
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